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Abstract 

In this study, we used a systems biology approach to Investigate changes In the proteome and metabolonne of shrimp 
hemocytes infected by the Invertebrate virus WSSV (white spot syndrome virus) at the viral genome replication stage 
(12 hpl) and the late stage (24 hpl). At 12 hpi, but not at 24 hpi, there was significant up-regulation of the markers of 
several metabolic pathways associated with the vertebrate Warburg effect (or aerobic glycolysis), Including glycolysis, the 
pentose phosphate pathway, nucleotide biosynthesis, glutamlnolysis and amino acid biosynthesis. We show that the PI3K- 
Akt-mTOR pathway was of central importance in triggering this WSSV-induced Warburg effect. Although dsRNA silencing of 
the mTORCI activator Rheb had only a relatively minor impact on WSSV replication, in vivo chemical inhibition of Akt, 
mTORCl and mT0RC2 suppressed the WSSV-induced Warburg effect and reduced both WSSV gene expression and viral 
genome replication. When the Warburg effect was suppressed by pretreatment with the mTOR inhibitor Torin 1, even the 
subsequent up-regulation of the TCA cycle was Insufficient to satisfy the virus's requirements for energy and 
macromolecular precursors. The WSSV-induced Warburg effect therefore appears to be essential for successful viral 
replication. 
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Introduction 

The Warburg effect, which was first described by Warburg iir 
the 1930s, is a metabolic rerouting used by tumor cells and cancer 
cells to support their high energy requirements and high rates of 
macromolecular synthesis [1,2]. In cancer cells, the main hallmark 
of the Warburg effect is aerobic glycolysis, in which glucose 
consumption and lactate production are both increased even in 
the presence of oxygen [3] . Several other metabolic pathways are 
also enhanced, including the pentose phosphate pathway (PPP), 
amino acid metabolism and lipid homeostasis. The Warburg effect 
can also be induced in vitro by some vertebrate viruses, including 
human papillomavirus (HPV) [4]; human cytomegalovirus 
(HCMV) [5,6], Kaposi's sarcoma herpesvirus (KSHV) [7] and 
hepatitis C virus (HCV) [8], and recently we reported an in vim 
Warburg-like effect that was induced in shrimp hemocytes by the 
white spot syndrome virus (WSSV; genus Whispovirus, family 
Mmaviridai) [9]. 



WSSV is a large unique, complex, dsDNA virus, and in shrimp 
hemocytes, its complete in vivo replication cycle takes 22-24 h 
[9,10]. Although over 90% of WSSV viral genes show no sequence 
homology to any other known genes, some of its genes are known 
to express at different times in its replication cycle, including the 
immediate early gene iel, the early gene DNA polymerase (dnapotj, 
the late structural protein gene vp28 and the very late DNA mimic 
protein gene icpU. Our previous study showed that WSSV 
induced the hallmarks of metabolic changes associated with the 
mammalian Warburg effect at the beginning of its genome 
rephcation stage (12 hours post injection [hpi]) [9]. However, since 
this was the first time that an invertebrate virus had been shown to 
produce this kind of effect, in the present paper, we look more 
closely at the global metabolomic and proteomic changes induced 
by WSSV in order to confirm that all of the interrelated metabolic 
effects seen in vertebrate cells are also found in the invertebrate 
Warburg effect. For the metabolomic study, we used hquid 
chromatography-electrospray ionizatioii-tandem mass spectrome- 
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Author Summary 

The Warburg effect (or aerobic glycolysis) is a metabolic 
shift that was first found in cancer cells, but has also 
recently been discovered In vertebrate cells infected by 
viruses. The Warburg effect facilitates the production of 
more energy and building blocks to meet the enormous 
biosynthetic requirements of cancerous and virus-infected 
cells. To date, all of our knowledge of the Warburg effect 
comes from vertebrate cell systems and our previous 
paper was the first to suggest that the Warburg effect may 
also occur in invertebrates. Here, we use a state-of-the-art 
systems biology approach to show the global metabo- 
lomic and proteomic changes that are triggered in shrimp 
hemocytes by a shrimp virus, white spot syndrome virus 
(WSSV). We characterize several critical metabolic proper- 
ties of the invertebrate Warburg effect and show that they 
are similar to the vertebrate Warburg effect. WSSV triggers 
aerobic glycolysis via the PI3K-Akt-mT0R pathway, and 
during the WSSV genome replication stages, we show that 
the Warburg effect is essential for the virus, because even 
when the TCA cycle is boosted in mTOR-inactivated 
shrimp, this fails to provide enough energy and materials 
for successful viral replication. Our study provides new 
insights into the rerouting of the host metabolome that is 
triggered by an invertebrate virus. 



try (LC-ESI-MS/MS) to identify and measure the levels of 
intracellular metabolites in shrimp hemocytes at 12 and 24 hpi, 
i.e. at the beginning and end of WSSV's genome replication cycle. 
For the proteomic profiling, we used label-free proteomics at the 
same time points. Since activation of the PI3K-Akt-mTOR 
signaling pathway is used by cancer cells and viruses to trigger 
the Warburg effect [1 1-13], and it is also required for the effective 
replication of vertebrate viruses [14—17], in the second part of this 
study, we use in vivo drug treatments to investigate whether WSSV 
also uses this signal pathway to trigger the Warburg effect. 

Results 

Global proteomic analysis of shrimp hemocytes during 
acute WSSV infection 

To understand the global changes triggered by WSSV infection, 
hemocytes were collected from PBS- and WSSV-injected shrimp 
at the genome replication stage (12 hpi) and the late stage (24 hpi) 
of the first WSSV replication cycle [9]. Using a label-free 
proteomic approach, 868 proteins were identified and quantified. 
Using a hierarchical clustering algorithm that grouped the shrimp 
samples by their protein abundance (Fig. SI), we found that 
WSSV-infected shrimp hemocytes had difiFerent proteomic 
expression patterns at 12 hpi and 24 hpi compared to the 
corresponding shrimp hemocytes collected from PBS-injected 
shrimp (Fig. SIA & SIB). No such proteomic clusters were formed 
by the hemocyte samples collected from PBS-injected shrimp at 
different time points (Fig. SIC), while two main clusters were 
formed by the WSSV 12 hpi and WSSV 24 hpi groups (Fig. SID). 
Two of the samples, 12-WSSV#1 and 24-WSSV#2, were not 
assigned to the corresponding cluster, and we therefore excluded 
these two mis-assigned samples from our subsequent analysis. (We 
note, however, that even when these two samples are included, the 
overall protein changes are only very slightiy different. Please see 
Table SI to compare the results obtained with and without the 
inclusion of these two anomalous samples). 



Global metabolomic analysis of shrimp hemocytes 
during acute WSSV infection 

To further understand the cellular responses after WSSV 
infection, we also used a global metabolomic platform to measure 
the metabolic changes in shrimp during WSSV infection. In this 
study, LC-ESI-MS data on over 100 metabolites were collected at 
12 and 24 h after WSSV- or PBS-injection. However, since we 
were interested primarily in host processes that are involved in the 
Warburg effect, we focused particularly on a limited number of 
important host pathways, including glycolysis, the PPP, nucleotide 
metabolism and the TCA cycle. Our metabolomic and proteomic 
data are given in Supplementary Tables SI and S2. Changes in 
these pathways at 12 and 24 hpi are shown in Figure 1 and are 
described in more detail below. 

WSSV infection enhances glycolysis at the WSSV genome 
replication stage (12 hpi) 

In mammalian cells, the two main pathways of carbon 
metabolism, glycolysis and the TCA cycle, oxidize hexose sugars 
to form ATP and NADPH, or else convert the same sugars to 
precursors of nucleotides, amino acids, and lipids. In shrimp 
hemocytes, WSSV infection at 12 hpi has previously been shown 
to increase glucose consumption and lactate production in ways 
that resemble the Warburg effect, but details of the intracellular 
changes in the carbon metabolism have not yet been investigated. 

In WSSV-injected shrimp hemocytes at 12 hpi, there was a 
significant increase (p<0.05) in the glycolytic pathway metabolites 
glucose, dihydroxyacetone phosphate (DHAP), glyceraldehyde-3- 
P, 3-phosphoglycerate (3-PG), 2-phosphoglycerate (2-PG), phos- 
phoenolpyruvate (PEP) and pyruvate (Fig. 1). There was also a 
corresponding increase in the protein levels of several glycolytic 
enzymes (Fig. 1). Despite the increase in glucose uptake, there was 
no evidence of lactate accumulation at the intracellular level. 
However, LC-ESI-MS revealed that lactate levels in the hemo- 
lymph outside the cells were significantly elevated at this time 
(Fig. 2A). By contrast, at the late stage of the WSSV replication 
cycle (24 hpi), except for pyruvate, all of the above-mentioned 
glycolytic metabolites decreased relative to the PBS-injected 
control (Fig. 1). Except for hexokinase (HK) and alpha-enolase 
(ENOl), the associated glycolysis enzymes were also less strongly 
up-regulated at this time. Meanwhile, there was now a significant 
accumulation of lactic acid inside the cells, whereas levels in the 
hemolymph returned to normal (Fig. 2A). Taken together, these 
results indicate that glycolysis was up-regnlated at the WSSV 
genome replication stage (12 hpi), but not at the late stage of the 
WSSV replication cycle (24 hpi). 

Apparent re-routing of the glycolytic pathway during 
WSSV infection 

At 12 hpi, we noticed that, even though glycolysis was 
enhanced, three glycolytic metabolites, glucose-6-P, fructose-6-P 
and fructose- 1 ,6-BP, remained unchanged relative to the PBS- 
injected controls (Fig. 1). The reason for this appears to be that 
glucose-6-P is being rerouted from glycolysis to the first step of the 
PPP. Although this putative rerouting still needs to be confirmed, 
e.g. by carbon flux analysis, this would explain why glyceraldehyde 
3-P, which is one of the end products of the PPP, was also 
increased. As Fig. 1 shows, many of the metabolic intermediates in 
the PPP were strongly up-regulated at 12 hpi even though there 
was only a relatively weak increase in the protein levels of the 
corresponding enzymes. However, although the protein level of 
glucose 6-phosphate dehydrogenase (G6PDH) was not detected in 
the present study, our previous data shows that there was an 
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Figure 1. WSSV induces the Warburg effect In the cellular proteome and metabolome of shrimp hemocytes at the replication stage 
(12 hpi) but not at the late stage (24 hpi). Changes In the levels of enzymes and proteins (ellipses) and metabolites (rectangles) relative to PBS- 
injected controls are color-coded to represent up- (red) or down- (green) regulation. Yellow represents no change. Colorless boxes and ellipses 
indicate that no data was detected. Protein data were collected from 3-5 pooled samples of 5 shrimp using quantitative label-free proteomics and 
expressed on a logarithmic scale. Metabolomlc data were collected from 5-6 pooled samples of 10 shrimp using LC-ESl/MS. Numeric values for the 
proteomic and metabolomlc data are given In Tables SI and S2, respectively. 
dol:1 0.1 371/journal.ppat.1 0041 96.g001 
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Figure 2. (A) WSSV infection increased the plasma lactate production at the WSSV genome replication stage (1 2 hpi), rather than at 
the late stage (24 hpi). The levels of lactate production were profiled and measured by LC-ESI-MS. Each bar represents the mean ± SD from four or 
five independent samples, with each sample pooled from 10 shrimp. The asterisks indicates a statistically significant difference in lactate production 
(p<0.05). (B) WSSV-induced phosphorylation of 4E-BP1 was suppressed by treatment with the mTOR inhibitors Rapamycin and Torin 1. Two hours 
before WSSV injection, shrimp were treated with volume-matched solvent (PEG), Rapamycin (RAP; 0.02 |xg/g shrimp) or Torin 1 (TR 1, 25 |ig/g 
shrimp). At 24 h post WSSV injection, twelve shrimp were selected from each group and divided into four sets (A-D, E-H, l-L, M-P). From each set, 
four pooled samples (3 shrimp in each pool) were prepared by collecting gill samples and extracting the total proteins. Each pooled sample was then 
subjected to Western blotting with antibodies to phosphorylated 4E-BP1-P^^''''* and actin. (C) Shrimp LvRheb is induced after WSSV infection. At the 
indicated time points after WSSV injection, total protein was extracted from the gills of individual shrimp and subjected to Western blotting. The 
expression of LvRheb was significantly induced at 24 hpi. ICPl 1, which is a major WSSV very late protein, was used as a proxy to indicate the WSSV 
infection state. Actin was used as an internal control. (D) Time series showing that LvRheb dsRNA treatment successfully silenced LvRheb expression 
in the hemocytes of WSSV-infected shrimp. Three days after LvRheb dsRNA injection, shrimp were injected with WSSV. Hemocytes were collected 
after the indicated number of days and subjected to cDNA synthesis and real-time PCR. Injections of enhanced green fluorescent protein-dsRNA 
(EGFP dsRNA) and phosphate-buffered saline (PBS) were used as dsRNA controls. Each bar represents the mean ± SD from four pooled samples with 
3 shrimp in each sample. An asterisk indicates a significant statistical difference between groups (p<0.05). (E) Gene silencing of shrimp LvRheb has no 
significant effect on the expression of the WSSV genes IE1 and VP28 during the first replication cycle (—24 hpi). The mRNA expression of IE1 gene and 
VP28 were used as proxies to indicate the WSSV infection state. Each bar represents the mean ± SD from four pooled samples (3 shrimp in each 
sample). An asterisk indicates a significant statistical difference between groups (p<0.05). (F) Gene silencing of shrimp LvRheb also has no significant 
effect on the number of WSSV genome copies during the first replication cycle (—24 hpi). Experimental conditions were as described above. The IQ 
Real WSSV Quantitative System was used to measure the number of copies of the WSSV genomic DNA. Each bar represents the mean ± SD from four 
pooled samples (3 shrimp in each sample). An asterisk indicates a significant statistical difference between groups {p<0.05). 
doi:l 0.1 371/journal.ppat.1 0041 96.g002 



increase in this enzyme's activity [9]. This suggests that the 
putative boost in the PPP at 1 2 hpi might be driven either by the 
increase in glycolysis or by changes in the activity of critical 
enzymes. However, we note that other explanations might also 
account for the observed accumulation of PPP metabolites, for 
instance, it is possible that the conversion of ribose-5-P to PRPP (5- 
phosphoribosyl-1 -pyrophosphate) might be blocked. 

At 12 hpi, we also observed an apparent boost in nucleotide 
metabolism, as suggested by increases in the levels of the purine 
biosynthetic intermediates in ATP synthesis (ADP, ATP, dATP) 
and dGTP synthesis (GDP, GTP and dGTP) (Fig. 1). Pyrimidine 
intermediates, including orotate as well as uridine, uracil, UDP, 
GDP and dTDP, were also up-regulated. This putative boost 
occurred despite a decrease in the levels of PRPP, presumably 
because this metabolite, which plays a number of roles in the 
biosynthesis of purine and pyrimidine nucleotides, was being rapidly 



consumed. As mentioned above, if conversion from ribose-.^-P was in 
fact blocked, this ^vould also contribute to the reduced lf;vels of PRPP. 

AU of the above changes would generate building blocks for 
macromolecular synthesis that could be used by the virus for 
genome replication. We further note that at 24 hpi, almost all of 
these increases in the PPP and nucleotide biosynthesis had 
dissipated. 

WSSV infection disrupts the Tricarboxylic Acid Cycle 

The end product of glycolysis, pyruvate, can be further 
converted into two metabolites, alanine and the important 
metabolite acetyl-GoA, which occupies a central position between 
glycolysis, the mitochondrial TGA cycle, beta-oxidation, amino 
acid biosynthesis and lipid biosynthesis. At 12 hpi, although 
pyruvate levels were significantly increased, surprisingly there was 
a decrease in the levels of acetyl-CoA (Fig. 1). The absence of any 
consistent increase in the levels of citrate, aconitate and isocitrate 
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suggests that the acetyl-CoA was being diverted (via citrate) from 
the TCA cycle into Upid biosynthesis. At 12 hpi, several other 
metabolites in the TCA cycle were significantly increased while 

their respective enzymes were decreased in parallel. A similar 
disruption was also observed at 24 hpi. Thus for instance, 
accumulation of the same four metabolites (i.e. succinate, 
fumarate, malate and oxaloacetate) was observed at both time 
points. Although these changes might have resulted from a 
blockage of this pathway, we also note that, at 12 hpi, even though 
glutamine levels showed no change, glutamate was significantiy 
up-regulated. This suggests an alternative explanation: i.e. that the 
second half of the TCA cycle is being driven by the conversion of 
glutamine to a-ketoglutarate via glutaminolysis. Input from 
glutaminolysis at 1 2 hpi would also explain why levels of succinate 
to oxaloacetate were increased even though the protein levels of 
their respective enzymes were all down-regulated. Lastly we note 
that both glutaminolysis and partial disruption of the TCA cycle 
are both well-known consequences of the Warburg effect. 

Enhanced production of amino acid metabolism during 
WSSV infection 

At 12 hpi, WSSV infection induced an increase in several of the 
detected amino acids, including alanine, histidine, tryptophan, 
glutamate, proline and aspartate (Fig. 1). At 24 hpi, amino acids 
were apparently still being synthesized, but the pattern had 
changed. For example, at 12 hpi, but not at 24 hpi, glutaminolysis 
seemed to drive the enhanced levels of glutamate, proline and 
aspartate. Meanwhile, although the high levels of aspartate at 
12 hpi were probably being used for de novo pyrimidine biosyn- 
thesis, by contrast, at 24 hpi, when the amounts of the nucleic 
acids had returned to control levels, the aspartate was evidendy 
being converted to asparagine and methionine instead. We also 
note that at 12 hpi, levels of pyruvate and alanine were both 
elevated, whereas at 24 hpi, even though pyruvate levels somehow 
remained high (despite the glycolytic pathway mostly shutting 
down), alanine levels had returned to normal. This appears to be 
due to the preferential conversion of pyruvate into lactate instead 
of alanine at this time. Taken together, these data suggest at 12 hpi 
WSSV infection strongly up-regulates amino acid metabolites that 
are useful for protein synthesis, and especially those that are 
involved in glutaminolysis and the synthesis of nucleic acids. This 
presumably benefits viral genome replication. We note however 
that this is only a tentative conclusion that still needs to be 
experimentally tested, and at this time we cannot rule out the 
possibiUty that the observed increases in these amino acids are 
simply due to down regulation or reduced activity of their 
corresponding enzymes. 

WSSV-lnduced activation of mlORCI is not specifically 
Rheb-dependent 

Recent studies have shown that numerous cancer/tumor cells 
depend on the mTOR signaling pathway to trigger the Warburg 
effect for efEcient cellular proliferation [11-13]. Although mTOR 
is found in two complexes, mTORCl and mTORC2, only the 
former is thought to be important for cell growth, proliferation and 
cellular metabolism [18,19]. Activation of the mTORCl pathway 
causes phosphorylation of its downstream target 4E-BP1, so 
phosphorylated 4E-BP1 is often used as an indicator of mTORCl 
activity [18]. 

Here we found that the level of phosphorylated 4E-BP1 in 
pooled samples of shrimp gills was elevated after WSSV injection 
compared to the PBS controls (Fig. 2B, lanes A-H). Further, when 
shrimp were injected with the mTORCl inhibitor Rapamycin 



(RAP) or the mTORCl /C2 inhibitor Torin 1 (TRl) 2 h before 
WSSV injection, the WSSV-induced phosphorylation of 4E-BP1 
was suppressed in all but one of the pooled samples (Fig. 2B, lanes 
I-P). From these data, we conclude that WSSV infection results in 
mTORCl activation. 

Our proteomic data also suggests that the mTORC 1 pathway is 
activated, as shown by the increased expression of several proteins 
in the mTOR pathway at 1 2 hpi, including the mTOR activator 
Rheb (Fig. S2A). We used Western blotting to confirm that Rheb 
expression was increased in WSSV-infec:ted shrimp at 12 and 
24 hpi compared to the control samples at 0 hpi (Fig. 2C). 

To further determine if the increase in Rheb expression was 
important for WSSV replication, shrimp were treated with Rheb 
dsRNA to silence Rheb expression before being injected with 
WSSV (Fig. 2D). Although all three groups showed a dramatic 
increase in Rheb niRNA expression from 0 to 12 hpi, Rheb 
mRNA expression was still much lower through to 48 hpi in the 
Rheb dsRNA-treated group compared to the PBS and EGFP 
dsRNA control groups. However, despite the successful silencing 
of Rheb expression, during the first WSSV replication cycle (0- 
24 hpi), the lEl and VP28 expression levels were only modestly 
inhibited in the Rheb dsRNA-treated group compared to the 
EGFP dsRNA-treated group (Fig. 2E), and there was only a 
correspondingly small decrease in the viral copy number (Fig. 2F). 
From these results, we concluded that the WSSV-induced 
activation of the mTOR pathway is pardy, but not completely, 
dependent on Rheb. This conclusion was further supported by 
Western blots that detected phosphor)4ated 4E-BP1 after Rheb 
dsRNA treatment (Fig. S2B), thus showing that mTORCl 
activation could still occur even when Rheb was knocked down. 

mTOR activation is critical for the WSSV-induced 
Warburg-like effect 

To investigate whether activation of the mTOR pathway might 
regulate the WSSV-induced Warburg-like effect, we pre-treated 
shrimp with Rapamycin and Torin 1 to respectively suppress 
mTORCl- and mTORC 1/C2- activation before injection with 
WSSV. At the WSSV genome replication stage (12 hpi), although 
lactate levels were elevated in the hemolymph of the non- 
suppressed WSSV-infected shrimp controls, there was no signif- 
icant lactate accumulation in WSSV-infected shrimp when 
mTOR activation was suppressed either by Rapamycin or by 
Torin 1 (Fig. 3A). Taken together, these data suggest that the 
mTOR pathway plays an essential role in triggering the Warburg 
effect. 

In shrimp hemocytes, the mTORCl pathway is 
predominant for WSSV gene expression, while inhibition 
of mTORCl alone is not sufficient to reduce WSSV 
replication 

Having shown that mTOR activation plays an essential role in 
triggering the Warburg-like effect during WSSV infection, we next 
explored what happens to WSSV viral gene expression and viral 
DNA genome replication when the Warburg-hke effect is 
suppressed by inhibition of the mTOR pathway. As shown in 
Fig. 3B, both Rapamycin and Torin 1 significantly inhibited the 
expression of WSSV viral genes in every stage of infection, from 
lEl, through DNA polymerase (DNA pol), to VP28 and ICPU. 
Inhibition of the mTORCl pathway alone was sufficient to 
account for most of the observed reductions, but the reductions in 
gene expression were further augmented when the mTORC2 
pathway was also suppressed. Although mTOR is known to be 
important for a range of cellular activities, including the activation 
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Figure 3. At 12 hpi, the WSSV-induced Warburg effect was suppressed by inhibiting mTOR, PI3K and Al<t. (A) At 12 h post WSSV 
injection, the accumulation of plasma lactate that was seen in the PEG-injected control group was suppressed by both Rapamycin (RAP) and Torin 1 
(TR1). Each bar represents the mean ± SD from four pooled samples of hemolymph (3 shrimp in each sample). An asterisk indicates a significant 
statistical difference between groups (p<0.05). The effect of Rapamycin and Torini on (B) the expression of the WSSV genes lEl, DNA pol, VP28 and 
ICP11, and (C) WSSV genome copy number. Data were based on pooled samples of hemocytes (gene expression) or pleopods (genome copy 
number), with all samples being taken from the same sets of shrimp as those used in Figure 2B. Specifically, the PEG results were from sets E-H, the 
RAP results from l-L, and the TR1 results from N-P. (Please note that data from the set IVI was excluded because phosphorylation of 4E-BP1 was not 
successfully inhibited in this set.). (D) Each bar represents the mean ± SD of the relative glucose or lactate concentration of five pooled samples of 
hemolymph (3 shrimp in each pooled sample) under the indicated conditions. Statistically significant differences are shown by 1-3 asterisks, which 
respectively indicate p<0.05, p<0.005 and p<0.0005. 
doi:1 0.1 371 /journal.ppat.1 0041 96.g003 



of several transcription factors [18,19], our results are also 
consistent with the following two ideas: First, that triggering the 
Warburg-like effect is essential for the successful expression of the 
WSSV genes and second, that although the niTORC 1 and 
mTORC2 pathways are both involved in triggering the Warburg- 
like effect, the mTORCl pathway is predominant. Meanwhile, 
synthesis of the viral DNA genome was only significantly 
suppressed by Torin 1 (Fig. 3C). One possible reason for this 
unexpected result is that blockage of the mTORC 1 pathway by 
Rapamycin does not completely shut down the Warburg effect 
because it still allows some amount of "leakage" through the 
mTORC 2 pathway. This leakage would also explain the observed 
significant differences between Rapamycin versus Torin 1 
treatment in three of the WSSV gene expressions (lEl, VP28 
and ICPll). 

Activation of the upstream PI3K-Akt pathway is also 
critical for the WSSV-induced Warburg-like effect 

Since mTOR activation can be stimulated by the upstream 
PI3K-Akt pathway [18], we next investigated the effect of 
inhibiting Akt activation. At 12 hpi, the Warburg effect, which is 
defined as an increase in both glucose consumption and lactate 



production, was only seen after pretreatment with PBS (Fig. 3D), 
(Taken together Fig, 1 and Fig, 3D show that there was an 
increased uptake of glucose from the hemolymph into the 
hemocytes,) This WSSV-induced Warburg effect was blocked by 
pretreatment with MK2206 (an Akt inhibitor). The effect was also 
blocked by LY294002 (a PI3K inhibitor that also inhibits mTOR) 
and by Torin 1 , There was no evidence of a Warburg effect at 
24 hpi under any treatment condition (Fig, 3D), In a separate 
experiment we also found that pretreatment with LY294002 
suppressed phosphorylation of 4E-BP1 at 24 h after WSSV 
infection (Fig, S2C), 

WSSV replication depends on the PBK-Akt-mlOR pathway 

To determine if the PI3K-Akt-mTOR pathway is critical for the 
promotion of WSSV gene expression and genome replication, we 
investigated the effect of the three inhibitors, LY294002, MK2206, 
and Torin 1, As shown in Fig. 4, treatment with 0,625-25 \lg of 
LY294002 per g of shrimp body weight led to a significant 
decrease in WSSV gene expression (Fig, 4A-C) and WSSV DNA 
genome replication (Fig, 4D), A second experiment with another 
batch of shrimp further showed that treatment with LY294002, 
MK2206 and Torin 1 led to significant and equal reductions in the 
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Figure 4. WSSV replication was suppressed wKien tlie PI3K-AI<t-mTOR pathway was inKiibited by LY294002, IVIK2206 and Torin 1. 

Shrimps were pretreated with LY294002 and samples were collected at 24 h post WSSV Injection. Graphs (A-D) show the expression of WSSV (A) lEl, 
(B) VP28, and (C) ICP1 1 in shrimp hemocytes and (D) the viral copy number in pleopods from the same groups of LY294002-pretreated shrimp (n = 6- 
1 0 In each group). Graph (E) shows the mean WSSV copy number in five pooled samples of pleopods (3 shrimp in each pooled sample) collected from 
another batch of experimentally infected shrimp after pretreatment with the indicated inhibitor. Bars labeled with different letters indicate 
significantly different values (p<0.05). 
doi:1 0.1 371 /journal.ppat.1 0041 96.g004 



WSSV genome copy number (Fig. 4E). These results suggest that 
WSSV infection may trigger the Warburg-like efTect by activating 
tiie PI3K-Akt-mTOR patiiway. 

To investigate the role of PI3K only, we used the selective pan- 
class I PI3K inhibitor BKM120. Although we were only able to 



conduct a pilot study, our results show that pretreatment with 
0.625 ng BKM120 per g shrimp significantiy reduced WSSV copy 
number (Fig. S2D). This provides further evidence that activation 
of the PI3K-Akt-mTOR pathway is important for WSSV 
replication. 
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Inhibition of mTOR activation by Torin 1 prevented the 
WSSV-induced Warburg effect in shrimp hemocytes 

To investigate whether the decrease in WSSV gene expression 
and viral genome copies in Torin 1 -pretreated shrimp (Figs. 3 & 4) 
might be due to cancellation of the Warburg-like metabolic 
changes, we compared the metabolic profiles of WSSV-infected 
hemocytes at 12 hpi either with or without Torin 1 treatment. 
Figure 5 shows that when shrimp were pretreated with Torin 1, 
the WSSV-induced Warburg-like effect was no longer observed at 
12 hpi (cf Fig. 1 and Fig. 5). In particular, in the glycolytic 
pathway, glucose was no longer accumulated, while the accumu- 
lation of glucose-6-P suggests that there was no rerouting into the 
PPP. Unlike Figure 1, Figure 5 shows no universal increase in 
nucleic acid synthesis, and this is also consistent with the observed 
accumulation of PRPP. Intracellular levels of lactate also 
increased. Meanwhile although the TCA cycle and amino acid 
biosynthesis were both up-regulated, the decrease in WSSV gene 
expression (Fig. 3) suggests that these changes alone are insufficient 
for WSSV pathogenesis. In the Torin 1 -pretreated shrimp, there 
was also no evidence of any WSSV-induced Warburg-like 
metabolic changes at 24 hpi (Fig. S3A). In control experiments 
where Torin 1 -treated shrimps were injected with PBS, at 12- 



24 hpi the metabolic intermediates were either down-regulated or 
remained unchanged (Fig. S3B, Table S2). 

Discussion 

Although only a limited number of vertebrate viruses (HCMV, 
HCV, Dengue virus, influenza virus, HSV and KSHV) have been 
studied using a systems biology approach [5,6,8,20-23], in the 
present study, several metabolic pathways that are often altered in 
cancer/ tumor cells and virus-infected vertebrate cells also seemed 
to be similarly affected by WSSV infection (Fig. 1 and Tables SI, 
S2). These altered metabohc pathways included glycolysis, the 
PPP, the biosynthesis of nucleic and amino acids, and the TCA 
cycle. 

The hallmark of the Warburg effect is aerobic glycolysis, i.e. a 
high rate of glycolysis accompanied by lactic acid production, 
despite readily available oxygen [3]. The effect is induced by 
several vertebrate oncogenic viruses in cultured mammalian cells 
[13], and in vertebrate cells infected with viruses such as HCMV 
[6]. Similarly, at 12 h post WSSV infection, i.e. the genome 
replication stage, glycolytic enzymes and intermediates were 
strikingly elevated and glucose uptake and lactate excretion were 
both increased (Figs. 1, 2A and 3D). In vertebrates, glucose 



Genome replication stage of WSSV infection (12 hpi) in shrimp pretreated with Torin 1 
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Figure 5. Pretreatment with Torin 1 inhibited the WSSV-induced Warburg effect at the genome replication stage of WSSV infection 
(12 hpi) in shrimp hemocytes. Two hours after treatment with Torin 1, shrimp were injected with PBS or a WSSV inoculum. At 12 hpi, 5 pooled 
hemocytes samples (10 shrimp per pool) were collected from each group. Cfianges in tfie metabolomic levels of the WSSV-infected samples relative 
to the PBS controls are color-coded as described in Figure 1. Figure S3 shows changes in the metabolome at 24 hpi. Numerical data for 12 hpi and 
24 hpi is given in Table S2. Changes in the metabolome for Torin-PBS versus PEG-PBS are shown in Fig. S3, with numerical data given in Table S2. 
doi:1 0.1 371/journal.ppat.10041 96.g005 
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transporter modulates this increased glucose uptake [24]; similarly, 
WSSV proteins interacted with glucose transporter during WSSV 
infection [25]. 

We previously reported that the activity of G6PDH was 
increased at the WSSV genome replication stage [9]. In the 
present study hexokinase (HK) levels were also increased (Fig. 1); 
therefore, we inferred that glucose was being rapidly converted 
into glycolytic intermediates. However, in the absence of an effect 
on the levels of glucose-6-P (its immediate downstream metabo- 
lite), we suggested that the glycolytic pathway was being rerouted 
into the PPP. Although we still lack carbon flux analysis data to 
show this directly, our present observations of significant increases 
in PPP intermediates and compounds used for nucleotide 
biosynthesis now provides further evidence that WSSV does in 
fact induced up-regulation of the PPP. 

In cancer cells, the PPP is usuadly up-regulated to balance 
ceUular redox conditions and to ensure an adequate supply of 
ribose-5-phosphate (R5P) for nucleotide synthesis, a rate-limiting 
step in cancer cell proliferation [26,27]. Similarly, viral replication 
requires biosynthetic precursors supplied by host cell metabolism, 
and synthesis of viral DNA creates a great demand for DNA 
precursors to be supplied through salvage reactions or de novo 
synthesis [22,28,29]. To ensure an adequate supply of nucleotides, 
viruses use several strategies, including expressing virus-encoded 
enzymes that boost nucleotide synthesis, and increasing the flux of 
the PPP [6,8,21,23,29,30]. With its large DNA genome 
(-307 kbp) and rapid replication cycle (22-24 h/cycle), WSSV 
appears to use both strategies. The WSSV genome encodes several 
viral enzymes related to nucleotide synthesis, e.g., dUTPase, 
thymidine kinase-thymidylate kinase (TK-TMK), and the ribonu- 
cleotide reductase subunits RRl and RR2, all of which are 
implicated in modulation of host nucleotide synthesis during 
WSSV infection [31]. We found that the PPP intermediates and 
the downstream nucleotides were up-regulated only at the WSSV 
genome replication stage (12 hpi) (Fig. 1). Even through WSSV 
induced up-regulation of the host's PPP enzymes at both 12 and 
24 hpi (Fig. 1), we reported that the activity of Glucose-6- 
phosphate dehydrogenase (G6PD), the rate-limiting enzyme of the 
pentose phosphate pathway, was increased at 12 hpi [9]. We 
therefore inferred that the putative increase in PPP activity was 
mainly being dri\'en not only by increased glycolysis, but also by 
enzyme activity during WSSV infection. Furthermore, we 
speculate that the observed down-regulation of PRPP (Fig. 1) 
occurred for different reasons at difierent times. At 12 hpi, the 
putatively increased throughput of both the PPP and nucleotide 
synthesis pathways would result in a high demand for PRPP. This 
would lead to low levels of PRPP, while conversely, high levels of 
this intermediate at 24 hpi could be explained by these two 
pathways shutting down. However, this explanation of the 
observed changes wfll need to be confirmed by fiirther studies. 

WSSV also induced an increase in the levels of the free amino 
acids, most of which are either derived from TCA cycle 
intermediates or converted from pyruvate and glutamine. 
Mazurek et al [32] reported that the increased production of 
amino acids during HPV infection resulted from both aerobic 
glycolysis and glutaminolysis. WSSV has already been shown to 
induce aerobic glycolysis, and Figure 1 now provides evidence that 
glutaminolysis is also induced: the levels of glutamate and its 
immediate downstream product, prohne, were both up-regulated 
at the WSSV genome replication stage, while glutaminolysis would 
also explain the observed replenishment of the partiaUy halted 
TCA cycle, and thus the increased levels of oxaloacetate and 
aspartate. Taken together, aU of these results suggest that WSSV 



infection regulates both of these pathways to enhance the 
production of amino acids. 

Acetyl-CoA, which can be converted from pyruvate, occupies a 
nodal position at the bifurcation of the anabolic and catabolic 
pathways [33], and it is involved in glycolysis, the TCA cycle, 
generation of nucleic acids and lipid metabolism. Biosynthesis in 
cancer ceUs consumes enormous amounts of energy and acetyl- 
CoA, resulting in a shortage of the latter for catabolic oxidative 
processes [33]. In cancc^r cells under the Warburg effect, shortages 
of acet)'l-CoA also result from a hypoxia-inducible factor (HIF)-1- 
mediated mechanism that prevents pyruvate dehydrogenase 
(PDH) from catalyzing the conversion of pyruvate to acetyl-CoA 
[2]. In contrast, Munger et al. [5] reported a dramatic increase in 
acetyl-CoA during the Warburg effect in HCMV-infected ceUs. 
However, in the present study, there was a significant decrease of 
acetyl-CoA at 12 hpi, but no increase in citric acid levels (Fig. 1). 
Therefore, we concluded that the WSSV-induced Warburg effect 
consistently resembled the Warburg effect in cancer cells, although 
there were some differences from the Warburg effect induced by 
other viruses. 

With such a marked decrease in the amount of acetyl-CoA 
entering the TCA cycle, it was surprising to see an accumulation of 
the TCA cycle intermediates succinate, fumarate, malate and 
oxaloacetate at 12 and 24 hpi. Two possible mechanisms might 
account for this. First, in proliferating glioblastoma cells, 
glutaminolysis, which ultimately results in the conversion of 
glutamine to lactate, provides an anaplerotic source for the TCA 
cycle and allows it to produce enough NADPH to support fatty 
acid synthesis [34]. Since glutaminolysis was also apparently up- 
regulated in the WSSV-infected hemocytes at 12 hpi (Fig. 1), this 
might likewise serve to replenish the TCA cycle and explain the 
observed accumulation of the downstream TCA cycle intermedi- 
ates. Another possible reason for the accumulation of these 
metabolites may be the down-regulation of related enzymes in the 
TCA cycle, such as succinyl-CoA synthetase (SCS), succinate 
dehydrogenase (SDH) and fumarate hydratase (FUMH) (Fig. 1). 
Recent studies suggest that SDH and FUMH behave as classic 
tumor suppressors [2]. When these enzymes are down-regulated, 
their respective substrates, i.e. succinate and fumarate, accumulate 
in the mitochondria and leak out into the cytosol, where their 
presence leads to an inhibition of prolyl hydroxylase enzymes 
(PHDs). Inhibition of the PHDs aUows the activation of HIF, 
which in turn causes pseudo-hypoxia and enhances both 
neovascularization and glycolysis within the cell [35]. Both 
glutaminolysis and enzyme reduction may be important at 
12 hpi, but as Figure 1 suggests, the downregulation of SCS, 
SDH and FUMH may be more important at 24 hpi l)C'cause 
glutaminolysis has apparendy returned to approximately baseline 
levels at this time. 

The phosphorylation of 4E-BP1 (Fig. 2B) shows that in the 
absence of the inhibitors Rapamycin and Torin 1, WSSV infection 
activates the P13K-Akt-mTOR pathway. Ivl(;amvhil(;, as shown by 
Figure 3A and Figure 5, pretreatment with the mTOR inhibitor 
Torin 1 disrupts aU of the above WSSV-induced metabofic 
changes. Taken together, these results suggest that the WSSV- 
induced Warburg effect is being driven by activation of the PI3K- 
Akt-mTOR pathway. This pathway is known to trigger the 
Warburg effect in cancer cells [11-13], and it is also used by 
human papillomavirus to trigger the Warburg effect and achieve 
successful replication [4,36]. 

When the Warburg effect was shut down by using PI3K-Akt- 
mTOR inhibitors, viral gene expression and viral copy number 
were significandy reduced (Figs. 3 & 4). We note that even though 
the TCA cycle was up-regulated in the Torin 1 pretreated shrimp 
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at 12 hpi (Fig. 5), the expression of viral genes and viral copy 
numbers was still repressed (Fig. 3). We infer from this that even 
when the normal source of energy (i.e. the TCA cycle) is activated, 
it still fails to provide sulTicient energy and materials to meet the 
demands of the virus, and that these requirements can only be 
satisfied by shifting the metabohsm into aerobic glycolysis. 

Our results show that when the PI3K-Akt-mTOR pathway is 
entirely blocked, i.e. by either an upstream inhibitor of PI3K or 
Akt (Fig. 4 and Fig. S2D) or by Torin 1 (Fig. 3), viral copy number 
is significantly reduced. We show here that the Warburg effect is 
also important for WSSV replication. However, we unexpectedly 
found that although Rapamycin was able to block the Warburg 
effect and the expression of viral genes (Fig. 3A, B), it did not have 
any significant impact on virus genome copy number (Fig. 3C). 
We hypothesize that these apparendy inconsistent observations 
may be due to the mTORC2 pathway, which is somehow able to 
provide an alternative route when mTORCl is blocked by 
Rapamycin. It would be interesting to see if a similar result could 
be produced by blocking mTORC2 only; unfortunately our 
attempts to achieve this with dsRNA silencing have so far been 
unsuccessful. 

Interestingly, silencing Rheb, which acts as a positive regulator 
of mTORCl signaling, does not significantly affect virus 
propagation during the first replication cycle (Fig. 2D-F). Two 
possible mechanisms might account for this. First, during WSSV 
infection, mTORC 1 activation might not depend exclusively on 
Rheb. Some vertebrate viruses, such as HCMV, can modulate the 
PI3K-Akt-mTOR signahng pathway via a surprisingly large 
number of alternative viral proteins [37-41]. It therefore seems 
possible that WSSV might likewise express one or more viral 
proteins that functionally replace the mTOR activator Rheb or 
otherwise activate the mTOR pathway. The second possibility is 
that mTORC2 can, to some extent, act to compensate for the 
suppression of mTORC 1 . As shown in Figure 3A, mTORC 1 is 
essential for inducing the Warburg effect, and it is also 
predominant in viral gene expression (Fig. 3B). Nevertheless, 
successful viral genome replication can be achieved by mTORC 2 
alone (Fig. 3C). Although these two mTOR complexes normally 
differ in their substrate specificity and function, during HCMV 
infection, mTORC 2 becomes involved in the phosphorylation of 
the mTORCl substrates 4E-BP1 and S6K [39]. Kudchodkar etal. 
[39] proposed that this specific alteration was due to structural 
modification caused by the addition of viral or cellular regulatory 
proteins to the mTOR complexes. If something similar is 
happening in WSSV, this suggests that there must stiU be other 
important WSSV proteins involved in the modulation of the 
PI3K-Akt-mTOR signaling pathway that remain to be discovered. 

In this study, in addition to providing a global metabolic view of 
the invertebrate Warburg effect during WSSV infection, we have 
developed a model of how WSSV may trigger these metabolic 
changes. A schematic is presented in Figure 6. In this model, the 
PI3K-Akt-mTOR pathway is activated after WSSV infection, and 
at the WSSV genome replication stage (12 hpi), the mTOR 
complexes mTORCl and mTORC2 are both involved. The 
Warburg effect and the anaplerotic effect of glutaminolysis both 
benefit the virus by increasing the availability of the building 
blocks and energy necessary to meet WSSV's requirements during 
genome replication. While it will be important to verify this model 
using techniques such as carbon flux analysis, which might be 
conveniently carried out in the relatively long-lived hematopoietic 
cells of an in vito crayfish culture [42], another interesting subject 
for future study will be to identify the putative viral factors that are 
involved in the activation of the PI3K-Akt-mTOR pathway. 



Materials and Methods 

Virus inoculum and experimental animals 

The virus used in this study was the WSSV Taiwan isolate. To 
prepare the WSSV stock, hemolymph was extracted from WSSV- 
infected moribund shrimp and subjected to centrifugation 
(10,000 xg). The supernatant was then diluted with phosphate- 
bufiered saline (PBS) (137 mM NaCl, 2.7 mM KCl, 10 mM 
Na2HP04, 2 mM KH2PO4), and stored at -80°C. The 
experimental inoculum was prepared from this stock by dilution 
(10 ) with PBS. Litopenaeus vannamei shrimp (mean wc'ight: 3—5 g) 
were obtained from the Aquatic Animal Center at National 
Taiwan Ocean University. Before the experiments, shrimp were 
cultured in a water tank system containing filtered seawater 
(30 ppm at 25 to 27°C) for 1-3 days. In tire WSSV challenge 
experiments, the shrimp were challenged with WSSV inoculum 
(100 |jl/ shrimp) by intramuscular injection. 

Ethic statement 

All of the shrimp used in this study were obtained from the 
Aquatic Animal Center at National Taiwan Ocean University. 
These animals were specifically raised for research purposes and 
Taiwan does not require any additional permit or permission. 
Because the experimental animals were invertebrates, no specific 
permits were required for this study, and there is no official 
recommendation for the use of shrimp for scientific purposes in 
Taiwan. Nevertheless all of our experimental procedures, includ- 
ing animal sacrifice, were designed to be as humane as possible, 
and all animals were treated so as to minimize suffering at all 
times. 

Antibodies 

Primary antibodies used in this study include phospho-4E-BPl 
(Thr37/46) (Cell Signaling; Catalog No. 2855), Rheb (Cell 
Signaling; Catalog No. 4935), and Actin (MUlipore). The antibody 
that recognizes the major WSSV late protein ICPl 1 was prepared 
in the lab as described previously [43]. 

Preparation of the mTOR inhibitors Rapamycin and 
Torin 1 

Rapamycin (sirolimus) stock was prepared by dissolving 
Rapamycin powder (Sigma-Aldrich Co.) in 99% ethyl alcohol. 
Before use, this stock was diluted with PEG solvent (0.25% 
polyethylene glycol, 0.25% Tween 20 and 0.15 M NaCl). Torin 1 
(Tocris Bioscience) was dissolved in dimethyl sulfoxide (DMSO) to 
provide a stock solution. Before use, this stock was diluted with 
PEG solvent. 

Effect of Rapamycin and Torin 1 on 4E-BP1 
phosphorylation, plasma lactate, viral genes and virus 
copy number 

To evaluate the involvement of the mTOR complexes during 
WSSV infection, shrimp were pretreated with Rapamycin 
(0.02 g shrimp) or Torin 1 (20 Hg/g shrimp) by intramuscular 
injection 2 h before being challenged. Control shrimps were 
injected with PEG only. At 12 and 24 h after the pretreated 
shrimps were injected with WSSV or PBS, four pooled samples of 
gills, hemolymph, hemocytes, and pleopods were collected from 
each group with each pooled sample being taken from 3 shrimp. 
Western blotting was used to measure the protein levels of 
phospho-4E-BPl in the giUs. Lactate levels in the hemolymph 
were measured as described below. Real-time PCR was used to 
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Figure 6. Schematic representation of the WSSV-induced Warburg effect and the involvement of the PI3K-Akt-mTOR pathway at 
the viral genome replication stage (12 hpi) of the first WSSV replication cycle. For the proteins Investigated In this study, red Indicates a 
positive Involvement and gray Indicates a partial Involvement. Blue Indicates Important upregulated proteins and Intermediates from previous 
studies: Glucose transporter (GLU1) Is from Huang ef al. [25], and G6PDH Is from Chen ef al. [9]. Elevated metabolic pathways are shown In pink. 
Pathway Inhibitors are Indicated In black boxes. Dashed lines Indicate Inferred regulatory mechanisms that were not Investigated In the present 
study. 
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measure virus gene mRNA expression and virus copy number. See 
below for details of these procedures. 

Effect of the PI3K-Akt-mT0R inhibitors LY294002, 
MK2206 and BKM120 

To evaluate the involvement of the PI3K-Akt-mTOR pathway 
during WSSV infection, shrimp were pretreated with inhibitors 
LY294002 (0.625-1.25 |xg LY294002/g shrimp), MK2206 
(0.625-1.25 ng MK2206/g shrimp) and BKM120 (buparlisib; 
0.15625-25 |ig/g shrimp) by intramuscular injection 2 h before 
being challenged. Stock solutions were prepared by dissolving 
LY294002 (Biovision), MK2206 (Biovision) and BKM120 (SeU- 
eckchem) in 10% DMSO, and these solutions were further diluted 
with PBS before use. Control shrimps were injected with 0.01% 
DMSO in PBS. At 24 h after the pretreated shrimps were injected 
with WSSV, hemocyte and pleopods samples were collected from 
6-10 individual shrimp in each group. The hemocyte samples 
were subjected to real-time PGR to measure the mRNA levels of 
WSSV lEl, VP28 and ICPl 1 as described below. Real-time PGR 
was also used to measure the viral copy number in the pleopods 
samples. 

Label-free quantitative proteomic analysis of the shrimp 
hemocyte by liquid chromatography tandem mass 
spectrometry (LC-MS/MS) 

At 12 and 24 h after shrimp were injected with PBS or WSSV, 
4—5 pooled hemocyte samples (5 shrimp in each sample) were 
collected from each group using an anticoagulant (450 mM NaGl, 



10 mM KGl, 10 mM EDTA, 10 mM Tris-HGl, pH 7.5). After 
centrifugation at 10,000 xg for 1 min followed by washing twice 
with Ix PBS, hemocytes were resuspended with 0.25 x PBS to 
extract total protein for LG-MS/MS based label-free quantitative 
proteomic analysis. The total protein content of the lysates was 
quantified using a Bradford protein assay kit (Bio-Rad) with bovine 
serum albumin (BSA) added as an internal quantitative standard 
for each analysis. The lyophilized hemocyte protein lysate (60 Hg) 
was resolubilized in an 8 M urea/25 mM ammonium bicarbonate 
buffer, incubated for 1 h at 37°G with 2 mM dithioerythreitol, and 
alkylated for 1 h using 2.5 mM iodoacetamide at room temper- 
ature. Each sample was then diluted with 25 mM ammonium 
bicarbonate to a fmal urea concentration of 1 M and Trypsin 
(Promega) was added (1:50 w/w). After overnight incubation, 
protease activity was quenched by acidification of the reaction 
mixture with formic acid solution (pH 1-2). Aliquots of the 
peptide mixtures were desalted and concentrated on a Cig- 
StageTip (Proxeon Biosystems), and then eluted with 50% 
acetonitrUe in 0. 1 % formic acid. 

The resulting peptide mixtures were analyzed by online 
nanoflow liquid chromatography tandem mass spectrometry 
(LG-MS/MS) on a nanoAcquity system (Waters, Milford, MA) 
connected to an LTQ^ Orbitrap Velos hybrid mass spectrometer 
(Thermo Fisher Scientific, Bremen, Germany) equipped with a 
PicoView nanospray interface (New Objective, Woburn, MA). 
After loading onto a 75-nmx250-mm nanoACQUITY UPLG 
BEH130 column packed with G18 resin (Waters, Milford USA), 
the peptide mixtures were separated at a flow rate of 300 nl/ min 
using a linear gradient from 5% to 40% solvent B (acetonitrile with 
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0.1% formic acid) for 90 min. The LTQ^ Orbitrap Velos 
instrument was operated in standard data-dependent acquisition 
mode, automatically switching between full-scan MS and CID 
(collision induced dissociation)-MS/MS acquisition. 

For the CID-MS/MS top20 method, full scan MS spectra (from 
m/z 350-1600) were acquired in the Orbitrap analyzer at a 
resolution of 60,000 (at 400 m/ z) and an AGC (automatic gain 
control) target value of 10^. The 20 most intense peptide ions with 
charge states S2 were sequentially isolated to a target value of 
5,000 and fragmented in the ion trap at 35% normalized collision 
energy, with an a( ti\ ati()n q of 0.25, 10 ms activation time, and 
minimum ion selection intensity of 500 counts. 

Progenesis LC-MS (Noiiliii(-ar Dynamics, version 3.0) was used 
for label-free quantification analysis. AH raw spectral files were first 
aligned and only high quality peptide features (charge state>l and 
isotopic pattern> = 3) were used for the extraction of ion intensity 
data from the MSI spectra. To adjust for system errors such as 
sample loading and intensity shift across LC-MS/MS runs, the 
peptide ion intensity was normalized using the Progenesis LC-MS 
robust mean, which was derived from the peptide log2 ratio 
distributions between a reference and the targeted LC-MS/MS 
run. The peak list generated from the qualified peptide features 
was used to search against a combined database that consisted of 
an in-house white shrimp database, the shrimp white spot 
syndrome virus database from NCBI, and the common Repository 
of Adventitious Proteins database downloaded from the Global 
Proteome Machine in the MASCOT 2.3 server (Matrix Science). 
To keep the false discovery rate below than 5 % (as estimated from 
the target-decoy database), only peptides with a Mascot ion score 
greater than 17 were included for subsequent protein analysis. 
Proteins were automatically assigned to functional group, and only 
these proteins that met both of the following criteria were 
reported: 1) the protein had the most peptide hits within its group, 
and 2) the protein included at least one unique, quantifiable 
peptide. Protein quantitation was based on the sum of the total ion 
intensity of the unique peptides. Lastly, to map the results into the 
biological network, MetaCore network software (GeneGO) was 
used for pathway analysis of the expressed proteins. 

Using liquid chromatography electrospray ionization 
mass spectrometry (LC-ESI-MS) to monitor the effect of 
Torin 1 on the metabolome of WSSV-infected shrimp 

In this experiment, 2 h before shrimp were challenged with 
WSSV or PBS, they were pretreated with PEG or Torin 1 (20 Hg/ 
g shrimp) by intramuscular inj('( tion to produce a total of four 
experimental groups: the PEG-PBS group, the PEG-WSSV group, 
the Torin 1-PBS group and the Torin 1-WSSV group. At 12 and 
24 hpi, 5-6 pooled hemocyte samples (10 shrimp in each sample) 
were collected from each group using anticoagulant as described 
above. After centrifugation at 800 xg for 1 min followed by washing 
twice with 1 x PBS, the hemocytes were resuspended with 0.33 x PBS 
and kept on ice for 10 min. The samples were then centrifugxxi at 
10,000 xg for 10 min, and 100% MeOH was added to the supernatant 
at a ratio 1: 2. After being centrifiiged again at 10,000 xg for 10 min, 
the supematants were lyophilized, dissolved in 35 nl ddH20 and 
subjected to LC-ESI/MS metabolomic analysis as follows: 

To enhance the detection of the carboxylic acid and organic 
phosphate signals, 5 |.ll aniline/HCl reaction buffer (0.3 M aniline 
[Sigma-Aldrich, USA] in 60 mM HCl) and 5 |J,1 of 20 mg/ml N - 
(3-dimethylaminopropyl)-N'-ethylcarbodiimide hydrochloride 
(EDC; Sigma-Aldrich, USA) were added to each sample of the 
hemocyte residue. Each mixture was vortexed and incubated at 
25°C for 2 h, after which the reaction was stopped by adding 5 jxl 
of 10% ammonium hydroxide. The aniline derivatized samples 



were then analyzed using an LC-ESI-MS system consisting of an 
ultra-performance liquid chromatography (UPLC) system (Ulti- 
mateSOOO RSLC, Dionex) and a quadrupole time-of-flight (TOF) 
mass spectrometer with an the electrospray ionization (ESI) source 
(maXis UHR-QToF system, Bruker Daltonics). The shrimp 
metabolites were separated by reversed-phase liquid chromatog- 
raphy (RPLC) on a BEH Cig column (2.1x100 mm, Walters). 
The LC parameters were as follows: autosampler temperature, 
4°C; injection volume, 10 (tl; and flow rate, 0.4 ml/min. After 
pre-starting with 1 % mobile phase B (0. 1 % formic acid in ACN) 
for 4 min, the elution started from 99 "'o mobile phase A (0.1% 
formic acid in ddH20) and 1 % mobile phase B (0. 1 % formic acid 
in ACN). After holding at 1% for 0.5 min and raising to 60% over 
5 min, mobile phase B was further raised to 90% in another 
0.5 min, held at 90% for 1.5 min, and then lowered back to 1% in 
0.5 min. The column was then equilibrated by pumping 99% B 
for 4 min. The acquisition parameters for LC-ESI-MS chromato- 
grams were as follows: dry gas temperature, 190°C; dry gas flow 
rate, 8 L/min; nebulizer gas, 1.4 bar and capillary voltage, 
3,500 V. Mass spectra were recorded from m/z 100-1000 in the 
negative ion mode. Data were acquired by HyStar and 
micrOTOF control software (Bruker Daltonics) and processed 
by DataAnalysis and TargetAnalysis software (Bruker Daltonics). 
Each metabolite was identified by matching with its theoretical m/ 
z value and with the isotope pattern derived from its chemical 
formula. The identified metabolites were quantified by summing the 
corresponding area of the extracted ion chromatogram, and 
metabolite sign;d levels were presented as the mean of the 5-6 pooled 
hemocyte samples from each experimental group at each time point. 

To investigate the WSSV-induced metabolic changes in shrimp 
hemocytes, the fold changes in the PEG-WSSV group were 
calculated relative to the PBS injection group (PEG-PBS group). 
To investigate the WSSV-induced metabolic changes in the 
mTOR-inactivated shrimp, the fold changes in the Torin 1-WSSV 
group were calculated relative to the PBS injection group (Torin 1- 
PBS group). Lastiy, the effect of Torin 1 pretreatment was shown 
by calculating the fold changes of the Torin 1-PBS group relative 
to the PEG pretreatment group (PEG-PBS group). Student's i-test 
was used to identify statistically significant changes. 

In vivo knocl<-down of LvRheb expression by dsRNA- 
mediated RNA interference 

Preparation of the dsRNA was done following Wang et al [44] . 
Briefly, first, the partial sequences (approximately 300-400 bp) of 
LvRheb and EGFP were genc'rated and amplifuxl l)y PGR for use 
as linearizing DNA templates. Next, the T7 promoter sequence 
was incorporated into these linearized DNA templates by using 
PGR with the following specific primer sets: Experimental group: 
LvRheb-dsT7F531/LvRheb-R882 and LvRheb-F531/LvRheb- 
dsT7R882; Control group: EGFP-dsT7F/EGFP-dsR and EGFP- 
dsT7R/EGFP-dsF (see Table S3 for details). The T7 RiboMAX 
Express large-scale RNA production system (Promega) was then 
used to synthesize the ssRNAs according to the manufacturer's 
instructions. The corresponding ssRNAs were mixed and annealed 
to become dsRNA by incubation at 70°C for 20 min, followed by 
slowly cooling to room temperature for 30 min. After purification 
and precipitation of the dsRNA by phenol/chloroform/isoamyl 
alcohol extraction, the dsRNA were quantified by UV spectro- 
photometer and verified by agarose gel electrophoresis. The final 
dsRNA products were stored at — 80°C before being used in the 
following in vivo experiments. 

For the gene silencing experiments, the experimental group was 
injected with LvRheb dsRNA (1 |J,g/g shrimp), while the control 
groups were injected with EGFP dsRNA or PBS only. To 
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determine the efficiency of the gene silencing for pooled hemocytes 
samples (3 shrimp in each pool sample) were collected from each 
group at the indicated time points. Total RNA was extracted from 
these samples, and cDNA was synthesized using Superscriptase II 
Reverse Transcriptase (Invitrogen) with Anchor-dTv primer 
(Table S3). Real-time PCRs were then performed to measure 
the expression levels of LvRheb and EFl-a with the following 
specific primer sets: LvRheb-qF/LvRheb-qR and EFl-a-qF/ 
EFl-a-qR. In uninfected shrimp, gene silencing was maximally 
effective at 3 days after dsRNA injection (data not shown). In 
subsequent experiments, the shrimp were therefore challenged at 3 
days post dsRNA injection. 

Quantification of WSSV gene (lEl and VP28) expression in 
LvRheb-l<nockdown slirimp 

For this knockdown experiment, shrimp were randomly di\ ided 
into 3 groups and injected with LvRheb dsRNA, EFGP dsRNA, 
or PBS. At 3 days post dsRNA injection, shrimp were then 
challenged with WSSV. Four pooled hemocyte samples were 
collected from each group at various time points (12, 24, 36, and 
48 hpi), with each pooled sample taken from 3 shrimp. Total 
cDNA was then prepared from each sample as described above. 
To quantify the relative expression of the WSSV iel and vp28 
genes, real-time PGR was performed with the specific primers lEl- 
qF/IEl-qR, VP28-qF/VP28-qR, and EFl-a-qF/EFl-a-qR (Ta- 
ble S3) using the Bio-Rad detection system with Brilliant SYBR 
Green QPCR master mix (Applied Biosystcms). Data values were 
calculated by the 2 ' mc-thod. Statistically significant differ- 
ences between groups were analyzed by Student's i-test. 

Quantification of the WSSV genome copy number in 
LvRheb-l<nockdown shrimp 

Four pleopod samples (3 shrimp in each sample) were also 
collected from of the above experimental groups at the same time 
points. The samples were subjected to genomic DNA extraction 
using a DTAB/CTAB DNA extraction kit (GeneReach Biotechnol- 
ogy Gorp.). WSSV genomic DNA copies were quantified using IQ, 
Real WSSV quantitative system (GeneReach Biotechnolog)-' Gorp.), 
which is a commercial real-time PGR based on the TaqMan assay. 

Determination of the concentration of hemolymph 
lactate in shrimp pretreated with Rapamycin, John 1, 
LY294002 and MK2206 after WSSV infection 

At 1 2 and 24 h post WSSV injection, 4—5 hemolymph samples 
(3 shrimp in each sample) were collected from groups of shrimp 
pretreated with LY294002, MK2206, Rapamycin, Torin 1 or 
PEG/PBS (control) without using anticoagulant. After being kept 
at 4°G for 12-16 hours, the samples were centrifuged at 13000xg 
for 15 min at 4°G, and the supernatants were transferred to new 
tubes. The concentration of glucose and lactate in the supernatants 
was then determined using enzymatic colorimetric test kits 
(Fortress Diagnostics Limited). 

Relative quantification of the expression of WSSV genes 
in shrimp pretreated with Rapamycin, Torin 1, LY294002, 
I\/1K2206 and BKM120 

After total hemocyte cDNA was prepared from all samples as 
described above, real-time PGR was performed with the specific 
primer sets lEl-qF/IEl-qR, DNApol-qF/DNApol-qR, VP28-qF/ 
VP28-qR, IGPl 1-qF/IGPl 1-qR and EFl-a-qF/EFl-a-qR (Table 
S3) using the Bio-Rad detection system with Brilliant SYBR Green 
Q_PGR master mix (Applied Biosystcms). Data values were 
calculated and presented as described above. Student's Hest was 



used to statistically analyze the Rapamycin and Torin 1 results. 
The LY294002 experiments used Tukey's multiple-comparison 
test (SPSS computer software) to evaluate statitiscally significant 
differences between experiemtnal groups. 

Quantification of WSSV genome copy number in shrimp 

pretreated with Rapamycin, Torin 1, LY294002 and MK2206 

Genomic DNA was extracted from pleopod samples, and the 
number of WSSV genomic DNA copies was quantified by the IQ^ 
Real WSSV quantitative system (GeneReach Biotechnology 
Gorp.) as described above. Data values were calculated, presented 
and statistically analyzed as described above. 

Protein extraction and western blot analysis 

Shrimp giU tissues were lysed in 0.33 x PBS with protein 
inhibitor and phosphatase inhibitor (Roche). Protein concentra- 
tions in each lysate were measured by Bio-Rad Protein Assay. 
Approximately 25 |J.g of protein lysate per sample were separated 
by \y/a sodium dodecyl sulfate-polyacrylamide gel electrophore- 
sis, transferred onto polyvinylidene fluoride (PDVF) membranes, 
blocked with 1-3% skim milk in Tris-buflered saUne with 0.1% 
Tween 20 (TBST) for 1 hour at room temperature, and then 
incubated overnight in primary antibody in TBST at 4°G. 
Following three extensive washes with TBST, membranes were 
incubated with horseradish peroxidase (HRP)-conjugated second- 
ary antibody (Santa Gruz) for 1 hour at room temperature. After 
three more washes with TBST, the signals were developed by EGL 
detection agents (Amersham) and detected using chemilumines- 
cence (Image Quant LAS 4000 mini). 

Supporting Information 

Figure SI Hierarchical K-means Clustering of ~800 hemocyte 
protein expression profiles obtained from a large-scale, high- 
throughput, label-free, quantitative LG-MS/MS analysis. (A), (B) 
At 12 hours and 24 hours post injection, the protein profiles of the 
PBS and WSSV groups formed two distinct clusters based on their 
log2 protein abundance, suggesting that the host cell protein 
pattern was markedly changed after WSSV infection at both time 
points. (C) Although there was no significant difference between 
the PBS groups at 12 and 24 hpi, the protein profiles of (D) the 
WSSV groups at 12 and 24 hpi formed two distinct clades, 
indicating that the host responses were different after WSSV 
infection at 12 and 24 hpi. Two of the samples, 12-WSSV#1 and 
24-WSSV#2, were not assigned to the corresponding cluster, and 
we therefore excluded these two mis-assigned samples from our 
subsequent analysis. 
(TIF) 

Figure S2 Proteomic data suggests that the mTOR pathway is 
activated at the replication stage (12 hpi) of WSSV infection. (A) 

Changes in the levels of enzymes and proteins (ellipses) relative to 
PBS-injected controls are color-coded to represent up- (red) or 
down- (green) regulation. Yellow represents no change. Colorless 
ellipses indicate that no data was detected. (B) WSSV-induced 
phosphorylation of 4E-BP1 was still detected even after Rheb was 
knocked down by Rheb dsRNA. Each lane shows the results for a 
pooled sample (n = 3) of total protein extracted from giUs and 
probes with antibodies against 4E-BP1-P^ IGPll and actin. 
(G) WSSV-induced phosphorylation of 4E-BP1 was suppressed by 
pretreatment with the inhibitor LY294002. Each lane shows the 
result for a pooled sample (n = 3) of total protein subjected to 
Western blotting with antibodies against 4E-BP1-P^^ and 
actin. (D) WSSV replication was significantiy reduced by 
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specifically suppressing using pretreatment with 0.625 [Ig/g 
shrimp of the selective pan-class I PI3K inhibitor BKM120 [45]. 
Data represent the mean ± SD of five pooled samples with each 
sample being taken from three different shrimp. 
(TIF) 

Figure S3 In Torin 1-pretreated shrimp, the Warburg effect was 
not seen either at 24 hpi in WSSV-infected shrimp or at 
12~24 hpi in PBS-injected shrimp. (A) Two hours after treatment 
with Torin 1, shrimp were injected with PBS or a WSSV 
inoculum. At 24 hpi, 6 pooled hemocytes samples (10 shrimp per 
pool) were collected from each group. Changes in the metabo- 
lomic levels of the WSSV-infected samples relative to the PBS 
controls are color-coded as described in Figure 1 . Numerical data 
for 24 hpi is given in Table S2. (B) Effect of Torin 1 pretreatment 
at 12 and 24 h post PBS injection. The metabolic intermediates in 
Torin 1-pretreated shrimps injected with PBS were either down- 
regulated or remained unchanged. Changes in the metabolome for 
Torin 1-PBS versus PEG-PBS at 12 hpi and 24 hpi are shown in 
color-coded boxes as described in Figure 1, with numerical data 
given in Table S2. 
(TIF) 

Table SI Global changes in the shrimp hemocyte proteome 

after WSSV infection. 

(DOCX) 
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